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ABSTRACT 



We study the competition between induced symmetry breaking (ISB), 
Cooper pairing (superconductivity) and chiral condensation at finite non- 
asymptotic density. Using a quantum mechanical model studied recently 
by Ilieva and Thirring, we analyze the expectation value of the fermion 
number density {a* a) in the presence of chemical potential and discuss rele- 
vance of the ISB phenomenon in the model. By Bogoliubov transformation, 
we obtain the ground state energy and find that a phase with both super- 
conducting and mean field gap is a true ground state. We also study the 
effective potentials of a two-dimensional field theory model for competitions 
between {'4>ip), {i/^ip) and (ip^t/j). We find that in the regime /i < ^c, where 
/ic is found to be around OAmp with mp being the renormalizaton point of 
chiral condensate, we have a chiral symmetry-breaking phase with almost 
zero fermion number density and in the case of ^ > //c the system sits in 
the mixed phase characterized by {ip^ip) ^ and {ipip) 7^ as long as the 
renormalization point, Mq, corresponding to a superconducting gap is not 
zero, Mo / 0. 



1 Introduction 



What happens to QCD at large density is a fascinating topic which is 
generating an intense activity in nuclear and particle communities. It seems 
that at an asymptotic density, the likely scenario is that diquarks condense 
into the matter as a consequence of which color superconductivity (CSC) 
may take place. For a recent review with references, see [||]. But it is not 
at all clear whether the density involved for CSC |^ is relevant to the dense 
matter that can be accessed by experiments and that we are interested in, 
namely neutron stars and heavy ion collisions that are currently studied or 
planned for the future. The question we would like to address here is: What 
happens to nuclear matter as density increases beyond the normal matter 
density toward the critical density for chiral phase transition? This question 
is not easy to answer since the coupling involved is not weak enough as what 
might happen at asymptotic density at which QCD weak-coupling property 
can be exploited but it is definitely more relevant to the physics we want to 
study in conjunction with experiments. We expect a priori that a variety 
of different phenomena will intervene in the regime of density that we are 
interested in and these may not render themselves to simple analyses. Some 
of them were discussed from a weak-coupling QCD point of view in Q. Here 
we address them from a different perspective. 

The issue in question was studied by Langfeld and Rho [Q] using a semi- 
realistic field theory model of the NJL type. Some of the results found there 
were novel and suggestive but certain restrictive features of the NJL-type 
models (i.e., no confinement) rendered them problematic. See for discus- 
sions on related matters. Here we would like to analyze the problem using a 
simplified solvable quantum mechanical model and a two-dimensional field 
theory model. Our study here is exploratory and we cannot make quanti- 
tative statements (for instance, our analysis is "blind" to color and flavor 
contents) as to what really must happen but our hope is to gain some quali- 
tative understanding of what might be taking place generically at a density 
greater than that of normal nuclear matter but much less than the asymp- 
totic density relevant to the color superconductivity. 

In Q, Langfeld and Rho suggested that when the isoscalar vector chan- 
nel (called the uj channel) becomes sufficiently attractive as expected in the 
presence of collective phonon modes of the pertinent quantum numbers as 
discussed in |p, then the baryon density can have a significant jump at some 
critical chemical potential in conjunction with generating low-mass excita- 
tions with the quantum numbers of the uj meson. Since this can be expressed 
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as a rapid increase - but not a discontinuous jump - in the ground state ex- 
pectation value of the time component of the oj field, the process was called 
"induced symmetry breaking" (ISB) associated with the broken isoscalar 
vector symmetry characterized by the "vaccum" expectation value of the 
time component of the vector current. We should stress that in medium, 
{ujq) is of course never zero, being proportional to density. Therefore {loq) 
cannot, strictly speaking, be used as a signal for a phase change. However it 
can have an anomalous increase at certain chemical potential reminiscent of 
an order parameter that can be associated with a phase transition. A sim- 
ple condensed matter model that illustrates this phenomenon is discussed 
by Langfeld 0. We shall use somewhat abusively this notion in this paper 
with the caveat in mind. 

A similar behavior of the baryon density vs. the chemical potential seen 
in the ISB has been also observed by Berges et. al |^ in a renormalization- 
group flow analysis of a linear sigma model in which it was found that 
quark number density shows non-analytic behavior at = 1.025Mq, with 
Mq ~ 2>lQ.2MeV an effective constituent quark mass, where the vacuum 
expectation value of the a-field, chiral condensate, vanishes. As discussed 
in ref.Q, this is easy to understand in terms of change in the number dis- 
tribution at that /i. The change in the number distribution can give a jump 
in the number density through /9 ~ / dk[e'x.]i (5{e{k) — fj.) + That there 
is a rapid increase in density at the transition point is intrinsic in the chiral 
phase transition. The point of ref.[^, however, was that the ISB is concur- 
rent with chiral symmetry restoration and that the presence of the ISB could 
postpone color superconductivity until the w-channel becomes ineffective at 
some high density due to weakened QCD coupling. At what density this can 
happen, we cannot say but the density regime involved could be relevant to 
the interior of compact stars. 

In this paper, we wish to address the problem in two aspects. First we 
take a quantum mechanical model studied by Ilieva and Thirring (IT) ^] 
and study, exploiting its solvability, whether and how the ISB manisfests it- 
self in the model. As stressed by IT, this analysis does not depend upon 
the dimension of the space considered. We shall discuss in section 2 the ISB 
phenomenon in the IT model and find that the model can have a phase in 
which both the ISB and fermion pair (Cooper) condensate co-exist and find 
that the ground state energy of the mixed phase is indeed lower than the 
normal state without such condensates. This model does not however ren- 
der itself to a discussion of what happens with chiral (fermion-antifermion) 
condensates. To investigate whether this kind of mixed phase is possible and 
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also how the chiral condensate figures in that phase, we resort to a soluble 
field theory model in two space-time dimensions studied by Chodos et al [^] 
and examine the effective potentials to see whether both condensates can 
coexist in the global minimum. The two-dimensional model is subject to 
the Mermin-Wagner-Coleman [O] no-go theorem for spontaneous symme- 
try breaking but we follow ref.[^] in considering this model in the large N 
sense a la Witten It would be possible to formulate the problem in ef- 
fective two dimensions dimensionally reduced from four dimensions, thereby 
avoiding the Mermin-Wagner-Coleman theorem. We find in section 3 that 
although gap equations exist in which both condensates are non-zero, the 
global minimum of the effective potential always occurs for the case when 
one or the other condensate vanishes as long as (V'V') and (V'^V') are con- 
cerned. On the other hand, we do find a stable mixed phase consisting of 
(il^ip) and (ip^il^) at some high density. Section 4 contains concluding re- 
marks. In Appendix we clarify the notion of grand canonical ensemble used 
in this paper and in pO[| . 



2 Ilieva-Thirring Model and the ISB 

Ilieva and Thirring |1C] studied the competition between {tp*ip) and (ipip) 
in a solvable quantum mechanical model (that we shall refer to as IT). Here 
we revisit their discussion in light of our objective as defined in the previous 
section. Since we shall use the notation and the results of IT, we first 
summarize them and then extract the information relevant to us. 

It was argued in [^] that two Hamiltonians are equivalent if they lead 
to the same time evolution of the local observables. This means that the 
effective Hamiltonian Hb 

Hb = J dp{a*{p)a{p)[u;ip) + Am{p)] + ^ABip)[a* {p)a* {-p) + a{-p)aip)]} (1) 
is equivalent to H = Hkin + Vb + Vm , 
H = Hkin + Vb + Vm 

where Vb and Va/ support respectively the pairing gap and the "mean- 
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field gap" Aa^ for k ^ cxo, provided tlie gap equations 

[a{k), Vb + Vm] = [a{k),HB - Hkin] (2) 

are satisfied. To solve the gap equations, IT choose the following potential 
describing an interaction concentrated about the Fermi surface 

VBMi^^P) = ^B,MS{k)S{p) 

with S{k) = L[Q{\k\-^ + e)-e{\k\-^-e)] (3) 

where Q{x) is the Heaviside function. With this potential and the additional 
assumption uj{p) = = ^ (with 2m = 1), the gap equations read 

W = Xb tanh ~ or = (5) 



w = ^[^i + ^M? + ^l (6) 

(where the pairing gap equation can have two solutions (|5|), the former non- 
trivial and the latter trivial) with the subsidiary conditions 

c\^i) + s\ii) = 1, 

2c{fi)s{fi) = AB{fi)/W (7) 

where c and s are real coefficients. It follows from (^ that sgnA^/ = sgnA^- 
It has been assumed that W > 0. From (|l|), we have the following thermal 
expectation value of the number density operator a*a, 

la-(p)a(p'))=S(p-p'){^-!^ + ^-pl>^}. (8) 

2.1 Ground-State Energy 

Let us now reformulate the IT model by performing Bogoliubov trans- 
formation explicitly and evaluate the ground-state energy which is related 
to the "effective energy" from which we can derive the same quantities. We 
write the reduced Hamitonian in the form 

H = Y^ ekaLflka + ^k^akolkO-iai + ^ V^Y 0^0^01,01. (9) 

kff kl kl 
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Introduce "order parameters" dk and mk by writing 

0-k«k = c^k + (a-kfflk - c?k) 
fflkflk = "T-k + (ak^k - "Ik) (10) 

and define 

Ak = -E^k^^i = -E^w(«-i«i) 
1 1 

Ak = -E^kr"^i = -E^kfK«i)- (11) 

1 1 

In terms of Ak and Ak, tlie reduced Hamiltonian (^) becomes 

= E ekaLak<7 - E^^k^ka-k + AkO-kOk - Ak^k) 

kcr k 

+ E(2Akakak- Ak^k). (12) 

k 

To diagonalize the Hamiltonian, we make use of the Bogohubov-Valatin 



canonical transformation |14| with real coefficients 

Ok = Ckfek + Sk^lk 

a-k = Ckb-k - Sk^k- (13) 
Substituting these new operators into we obtain 

H = «k[(ck - 4)(f'k^k + fe-k^-k) + 2ckSk&k^-k + 2ckSkfe-kfek + 2sl] 

k 

+ E[ckSk(Ak + A^)(6^6k + 6*_k&-k - 1) - (Ak4 - Al;4)&kfelk 

k 

+ (Ak4 - Aj;4)6_k&k + Ak4] + E[2^k(4f'k^k + CkSk^k^-k 

k 

+CkSkfo-k^k + s^fc-k^lk) - AkW-k]. (14) 

The diagonalization is effected by demanding that the coefficients of ^'k^lk 
and 6_k^k vanish 

2ckSkek + A^4 - Ak4 + 2AkCkSk = 0. (15) 
Multiplying AJ^/4 and defining x = ^AJ^, we get two solutions 



ek + Ak)± V(^k + Ak)2 + A2. (16) 



We take (+) sign here to get a stable minimum energy solution. Then 
assuming Ak to be a real quantity, we have two equations for Ck and Ski 

£k _ Ej, - (ek + Ak) 
Ck Ak 

4 + 4 = 1 (17) 



where £'k = y (ek + Ak)^ + A^. Solving the equations, we have 

4 = 5(1-2^). (18) 

To obtain the (coupled) gap equations for Ak and Ak, we rewrite Ak and 
Ak in terms of the new operators defined in (|T3|), i.e.. 



Ak = -Y.VMCiSi{l-blbi-bUb-i) 



1 



Ak = -E^kY(4 + 4&r&i-4^r&i) 
1 

= - II ^kl [^I3{E,\) + I + e-/3(Bi-M) + J ^^^^ 



1 



where we have used = l/(e'^(-^'~'^) + 1). To compare ([l9|) with 

the gap equations of IT model, we take = — 2AB5|k|,^'^|i|,^i = 
— 2AM^|k|,^^|i|,^ with 5afi being Kronecker delta, Ak = ^B{^J) and Ak = 
/S.M{^J) and get 

E{^l) = As tanh^^^^^^l^ or As (^) = 

Am(m) = 2XM[^f,^j,lJ;'^l^ ^ ^ + + i^- (20) 

We can easily read off the ground state energy U from (11^ 



U = E[2ek4-2AkCkSk + Ak4- Akmk + 2Ak4] 

k 

= ,-EM + A„M + ^-%M pi) 

IAb ^Am 
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where the summation has been performed with d\^^^ in the spirit of the 
IT modelQ. We are now in a position to analyze various cases of interest. 

2.2 Ab = 

This case corresponds to the simplest solution of gap equations for all 
values of Xm, Ab, fi. 



Case I-i : /u + Am > 
In this case, 

E- fi = Am 



2\m 



1 + e/3^M 



. f 2Xm for Am < , . 

^ { ^ ^ Q , (22) 

as /3 — > oo. Since the phase with Am = 2Am has a non-zero positive ground 
state energy U = — 2Am foi" Xm < 0, it does not interest us. We will not 
consider it anymore. Then, ^ becomes (after integrating over p') 

{a*{p)a{p)) = for Am = (23) 

where e = p^. Note that this phase is nothing but the normal state ( 
As = and Am = 0). 

Case J-ii : /i + Am < 
To satisfy the condition /x + Am < 0, Am should be negative and therefore 
Am < 0. Noting that (et + Ak)/£'k = — 1 when we take Ak(AB) = and 
Ck = in (H); we have c^(^) = and s^(^) = 1. Therefore we have 

E + ^i = -Am 

2Xm 



A 



M 



I + e/3(AAf+2At) 



. 2Am for Am + 2^ < , . 

^ forAM + 2^>0 ^ ^ 



^In the usual BCS theory, the summation is replaced by the integration, « ""^ 
with hujc, a typical phonon energy. But since we are considering the very special case of 
an interaction concentrated about the Fermi surface, we impose the Kronecker delta in 
the summation. 
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as (3 ^ oo. Noting that the second solution Am = cannot satisfy the 
condition + Am < 0, we have at a given energy e{p) = 

{a{p)a{p)) = , e + Am < 0, 

= i + e^U.-,) ^^ + ^M>0 (25) 

for Am = 2Am- The corresponding ground state energy U = 2[fi + Am) 
is negative as long as /u < |Am| and therefore this phase is energetically 



favorable compared with the normal state of matter. Since the Am in (25) 
corresponds to (Vq) in eq.(2) of ref.[^, one may be tempted to conclude that 
it is a signal of an ISB. But this can not be a candidate for an ISB because 
Am in the IT model cannot be associated with spontaneous breaking of a 
symmetry and furthermore Am is a constant independent of the chemical 
potential. 

2.3 Ab ^ 0, Ab < 

In this case, we cannot expect the Cooper pairs to condensate if Am = 
as one can see clearly from the gap equation (^). If we take Am = and 
therefore Am = 0, the gap equation becomes 



+ A| = Ab for Ab / and /3 ^ oo (26) 

which is contradictory with the conditions we are starting with. 

It follows that the parameters under consideration must respect the fol- 
lowing conditions |1C]: 

Ab < 0, E<^i, Am < 0, Am < 0. (27) 

Since we are interested in the physics of cold dense matter, we solve the gap 
equation taking /3 — > oo and get 0, 11] 



E = — Ab, Am = Am- ^ ^ 



Ab + Am 



Ab 



Ab = ±T — ^a/(Ab-^)(Ab + 2Am + m) . (28) 
Ab + Am * 

Note that we have solutions only when jj, > E = — Ab- To investigate the 
competition between the mixed phase given by (p8|) and the Aj\/ ^ phase 
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in /-ii, we take the following values of coupling constants with arbitrary 
dimension 



Si 


\b = 


-5, 


Am = 


-5 


S2 


Xb = 


-5, 


Am = 


-10 


S3 


Xb = 


-5, 


Am = 


-2 



(29) 

and find the following critical chemical potential fic 



Si 




= 5 


S2 




= 7.6 


S3 




= 2, f4 



5. (30) 

In the case of Si and 52, our system sits in the phase with = and 
Am 7^ 0, /-ii, in the regime fi < fic and rolls down to 7^ and Am 7^ 
above Hc- The case S3 is quite puzzling in the sense that in the region jjL < fi]. 
our system sits in the phase with = and Am 7^ 0, for /ij < /i < /x^, 
we have the normal phase (Ag = and Am = 0) and in the region fj, > fx^, 
our system sits in the A^ 7^ and Am 7^ phase. But in all three cases, 
we recover the normal phase at a high chemical potential. Si : fj, > 8.5, ^2 : 
fj. > 12.3, S3 : fi > 9.7. The phase diagram with the parameter set Si is 
shown in fig. [l|(a). To discuss a possible realization of ISB at the critical 
potential ^c, we note that in the mixed phase U{1) symmetry of the IT 
model is spontaneously broken to Z2 by the superconducting gap A^. Thus 
one may expect an ISB phenomenon to take place at the critical chemical 
potential but one should note that chemical potential does not break U{1) 
symmetry explicitly at the level of an action. One of the essential points of 
ISB is that a small explicit breaking of a symmetry by chemical potential or 
external magnetic field at the level of Lagranigan leads to drastic changes 
in the vacuum structure of the theory and thereby gives a jump in number 
density or magnetization [^. In addition, the mixed phase in the IT 
model does not possess the exclusive feature that the presence of the ISB 
expel color superconductivity found in ref. Q. Thus we are led to suggest 
that the ISB phenomena may not be a relevant notion in the IT model 0. 

2.4 Ab^O,Xb>0 



In this case, the solutions for A^ and Am are the same as those in (28) 



but E = Xb with a restriction Ab + /x + 2Am > in the limit /? — > cxd. We 
^We are grateful to K. Langfeld for emphasizing this point to us. 
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find from the ground state energy given by tlie parameter clioice ( |29| ) that 
as long as < //c = As, we have a mixed phase, 7^ and Am 7^ 0, 
and that in the regime /i > ^c, our system sits in the normal phase. The 
resulting phase diagram is depicted in fig. |l] (b). 



Ab=0 




Ab=0 






Am=0 



(a) 



(b) 



Ab=0 
Am=0 



Figure 1: The phase diagram calculated in IT model in cold dense matter. Here 
energy dimension of chemical potential is arbitrary. 



3 Field Theory Model For Competitions Between 

{ijj^p), {ipip) And {ip^ip) 

In this section, we shall study a field-theory model by generalizing the 
model considered recently by Chodos et al |]l5| to one that contains a number 
density field corresponding to (V'^V')- The (l+l)-dimensional toy model we 
shall consider is defined by the Lagrangian 

+2G2[Vi(^)75^(^')][V^^*^75V'^^'^] - 2F2[V;(^)7^V'^*^][V'^^'HmV'^^'^] (31) 

where i runs from 1 to N. The last term is the term we have added to 
the model of Chodos et al. We shall do a large approximation. We do 
a Hubbard-Stratanovich transformation using the auxiliary fields, (j), B and 
Vfj_ by adding the term 

+^K + F^^^^,i^^^f. (32) 
The resulting Lagrangian V = L -\- AL is 



10 



+Beo.pi>^J^^i^f - B^eo^pi^^^i^f ■ (33) 

We can now integrate out tjj and ip^ to obtain the effective action (modulo 
a constant) 

dM-^ - + p2^i) - 2Trln(ATA) 

-^TVln[l + M2(AT)-V2A-V2] (34) 

where we have assumed that B, B\ (j) and Vq are constant since we are inter- 
ested in an effective potential with translation invariance and have defined 

A = ido + ia^d^ — /x — 0(Ti 
= —ido — ia^ds — A — 4>(^i 

with fl = fi — 4Vb. 

Factoring out coming from the flavor trace, we write g'^N = A, G'^N = 
f and F'^N = 2a and define Veff as F^ff = -N{J d'^x)Veff: 



+ + ^^^^ 



(1) 



where v})j = i[Trln(ATA) + Trln(l + M2(AT)-1(72A-V2)]. Then the gap 
equations are obtained by 

from which we get 



2A J (27r)2 D 

1 _ dv}j} _. f d^k [k^ -kl-fl^-M^ + <^2] 



K dA'P J (27r)2 D 

1 _ dv}}l _ r d^k AHikj + kf-f,^-M^ + cp^] ^^^^ 



a dV^ J (27r)2 VqD 
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where D = [k^ - kj - jl'^ - M'^ + cf?]^ - Alcf^kl + /i^fcf - (t?kl] with ko 
ko + ie sgn^o- After the ko integral, we have 

_1_ ^ 1 ^^^^^ 1^1^ AP + fi^ ^ 1 -\_ 

- = ^ dki[jr + Tr+ I (7- -7-)] 

K 87r7_A k+ k^ yM202 + /i2(^2^02) /c+ /C- 



a Vb27r7-A k+ k^ ^ M^^ + fx^k^ + (f)"^) ^+ ^ 

where A is a cutoff to regularize logarithmic divergences 



and k± = \InP + 4)"^ + ft? + kl± 2y'M2,/.2 + fi'i{kl + 02). 

^6 



By integrating with respect to cfP' , M'^ and Vq, we have the unrenormal- 



ized V}j}, 

Kff = -^l dki[k+-k^ + c{ki)]+C{fi) (38) 

where c{ki) and C(/u), the value of which are determined later, are constants 
of integration. In the unrenormalized effective potential in free space, there 
can be quadratic and logarithmic divergences. Since free-space renormal- 
ization will eliminate all deivergences even in medium, we will cancel out 
the quadractic diverence by choosing a suitable value of c{ki) and define 
renormalized coupling constants to eliminate logarithmic divergences. To 
be explicit, we take /i = 0. Then the V^jj. is given by 



V}jj = dki M2 + 02 + 16^02 + ^2 + 2y/m202 + 16Vo\kf 



+ y M2 + 02 + 16^02 + kf- 2^M202 + IQV^ikf + 02) _ 2ki] (39) 

where we have chosen c{ki) = —2ki to cancel the quadratic divergence, see 
(1^ and C{fi = 0) = 0, see 



3.1 = 0, = 

Since it is not so easy to integrate (^) explicitly to separate the infinities 
from finite quantities, let us first confine ourselves to the case of = 0. 

^In Ref. it is shown that although a solution to the gap equations exists in which 
both condensates, (tpil^) and (i^tp), are non- vanishing, the global minimum of the effective 
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The y}ff is given by 



v}j}{Vo,M) = dkii^M^ + (41/0 + kif + \Im^ + (4^0 - - 2fci] 

= -^[(A + 4yo)^(A + 4^0)2 + M2 + (A - 4Fo)^(A - ^V^? + 



Air 



4Fo + -/lei^TM^ 



A - 4^0 + V(A - 4^o)^ + ^ 



+M^ln " ^ \ 2A^]. (40) 

-4^0 + JWV^ + M2 



The unrenormahzed effective potential is given by 

Veff{Vo,M) = M2(i - ^) - Fo'(i + ^) 
At 47r a TT 



1 ^^^^2 A + 4^0 + V(A + 4^0)^ + 
4^0 + ./i6K2 + M2 



^^,^^A-4W(^^M^]. (41) 
-4^0 + Vl6Fo^ + M2 



We define the renormahzed coupUngs kr and aj? as 



and get 



1 _ 1 + 8 

Oi? a TT 



1 11, A + Avo + JiA + Avof + M^ 
— = —[In , = 



A-4t;o + J(A-4t;o)2 + M2 

+ ln ^\ — -1]- (43) 

-4?;o + V16i;2 + 



potential always occurs for the case when one or the other condensate vanishes in free 
space {jj, = 0) except for one very special case which we will not consider here. If we take 
naively this result, setting = is not so bad an approximation. 
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The renormahzed effective potential is now found to be 

Vj}, = MH^-^)-V„^^ + ^l.^_. (44) 

''^ KR 2tt' ^ or Att ^ ^ 

Note that Vq and M are not coupled to each other. Solving the gap equations 
for Vq and M, we obtain 

Vo = 

= M^e^~^ or 0. (45) 
At the solutions, the effective potential becomes 
V^jf = for Vo = 0, = 

V3f = —r^e for Vo = 0, 0. (46) 

So the phase with Vq = 0, / is energetically favored and is a global 
minimum when or < 0. 

3.2 = 0, yU^O 

In this case, the renormahzed effective potential is given (up to a con- 
stant) by 

K^/ = -h- - + ^^Vo + ^ In |J + C(M) (47) 

kr In ur zvr vr 47r Mq 

where C(^) wih be fixed by the condition that V^jj{M = 0, = 0) = 0. 
We will make use of the condition V^^{M = 0, Vq = 0) = throughout this 
paper. Then we have the renormahzed effective potential, 

= M^(± - ±, - VSI- + + f in (48) 

kr 2tt ur vr 47r Mq 

As far as is concerned, our effective potential is the same as that in Ref. 



[15|. Note that since Vq and M do not couple, it is easy to find solutions of 



the gap equations. We find 



Vo = 2— /i, = or (49) 

TT 

Vo = 2—fi, = M^e~^. (50) 

TT 
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At the solutions, the effective potential becomes 



4 2 
Klf = -^aRfi"^ for Vo = an-ji, = 

K'lf = - ^e'"^ for Vo = aR-fi, / 0. (51) 



From (51), we can see that the phase with M 7^ is energetically favored 



regardless of the behavior of Vq. Thus we expect there could exist a phase 
where superconductivity and mean fields coexist in the case of Mq 7^ at 
finite density. Since we are expecting superconductivity at high density, it 
is plausible that we have Mq 7^ with Mq being small at some relevant 
density. 

3.3 M = 0,// = 

In this case, V^jj- takes the form 
1 /-A 



K// = -^J^ \/^^+^ + 4^0 I + I - m I -2h]. (52) 

We can consider two cases: 

• |<^| >4|yo|: 

In this case, the effective potential 1/ ^j. is independent of Vq. Using the 
exactly same method adopted in the previous section, we have 

1 11 a + ^/attW 1 

- = — + -ln (53) 

A Xr TT 00 



1 _ 1 

a ur 



(54) 



1 3.1 



where Xr is defined by q^I' |m=Mo,</>=0o~ X^- Since a and aR differ by 
a constant, we will use, hereafter, a instead of ur just for simplicity. The 
solutions of the gap equations are given by 

^0 = 0, </> = (56) 
or Vb = 0, ct) = (l)oe~^. (57) 
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It can be easily seen that the phase with Vq = and 7^ 0, at which 

2— — 

Vjjj = —-^4>Qe , is favored. This comes as no surprise. 

• m<MVo\: 
We find 

1 1 2 16|?;o| 



a an TT Jig^g _ ^2 



1 I I A+JA^ + ^I 1 <^2 
- + - In ^ , + " 



1 4|yo| + Vl6V^o' 



+ — 6^ In- 



2^ 4|^;o| + v/l6t;g-0g 



-_|yo|y^l6y,2_^2. (58) 



Solving the gap equations, we find that the only possible solutions are Vq = 
and (/) = where ^^ff = 0. So at zero chemical potential, the system is again 
characterized by Vq = and 7^ 0. 

3.4 M = 0,^^0 

For /i <| ^ I, the effective potential is equal to (55). In this case, our system 
is characterized by 



For /X >| (/> I, we have 



^0 = 0, 4> = 4>Qe ^n. (59) 



p J_ 2_ (go) 
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(Recall that ft = fi — 4Vb-) The possible solutions of the gap equations are 
calculated to be 



1 1-^ 



JL - A' 

2a ^ 

,2 



r = A^-i^^o- (61) 
At the solutions of the gap equations, we have 

^.Jz-jt-^H^J^W-^-^W. (62) 



<4|yo|: 



, 1 ,2i I A I - 0^ , 1 2 

2^ 16i;2-<^§ + 4K|^16i;2-<^2 27r 

First, we consider the case of 

<^ = and = (64) 

which is one of the solutions of the gap equations. In this case, the effective 
potential becomes 

Ve% = -^(16 + ^)Vi + ^Vofi. (65) 
For /i > 0, we have the following solution for the gap equation 

2a ^ 
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In the case of /i < 0, we have 

= -^'-^^yi (67) 

To see which phase is energeticahy favorable, we take (po = and Mq ~ 
0. 

In free space (/i = 0), we have two phases: 



y« = - rrijpe for = and (/> / 
V;^^ = --^e for = and M / 0. (68) 



An 

Since Mq will be zero or small, our system sits in the phase with Vq = and 
(j) ^ 0. In free space, therefore, we have chiral symmetry breaking but no 
sup er conductivity. 

At finite density, the situation becomes more complicated. We have six 
sets of solutions of the gap equation: (H), (|9|), (||), (g), (|6D and (|^. 



Now let us investigate the competition between the phases characterized 
by (|50|), ( ^9| ) and (|6^). The effective potentials given at their ground-state 
positions are 

a ' IT 

II : K// = ^ (ID 

m : = ^t^m' - (H). (69) 

a IT 

Comparing / and ///, we find that our system will be in the mixed phase 
(I) characterized by (ip'^ip) 7^ and (V'^) 7^ as long as Mq 7^ 0. From 
the condition that the phase / is the global minimum of the potential ([48D, 
^ + f < 0, we find — 0.1257r < a < 0. It is expected that the phase //, which 
is valid when /i < \(p\, becomes the absolute minimum of the potential at 
low density. To see the competition between I and // explicitly, we neglect 

^ Since it is not plausible to have superconductivity in free space, Mo should be zero 
or very small, if any. We shall assume, however, that we can have a small but non-zero 
value of Mo at finite density 
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a term with Mq in the phase I and take Xr = tt and a = — O.Itt. Then 
we find ~ 0.4. In the regime // < /ic, we have the chiral symmetry 
breaking phase (//) with almost zero fermion number density defined by 

Pf = III the case of /i > the system sits in the mixed phase (/) 

with fermion number density 



8aR , 1 1 

PF = o — z A* ^ith — = - + 

vr^ — 87ra/j an a vr 

0.25/i for a = -O.lvr. (70) 



The fermion number density is depicted in Fig. |2|. We should, however, 
consider competitions with other possible phases. To investigate numeri- 



cally the phases characterized by (plD, (66) and (p^), we take the following 
parameters 

\r = TT 

a = — vr, — O.lvr 
— = 10, 1.0, 0.3 . (71) 

With this parameter choice, we plot the renormalized effective potentials 
near the minimum-energy positions (i.e., solutions of the gap equations), 
(SI), ( |66[ ) and (|67|). We find that the phases characterized by (|6T1), ( |6^ ) and 
(67) correspond to the maximum or saddle point of the renormalized effective 
potential or does not satisfy some constraints, for example 4|Vo| > \4>\, near 
the minimum-energy positions and therefore those phases are unstable or 
an unphysical "vacuum." 



4 Conclusion 

We studied the competitions between induced symmetry breaking (ISB), 
Cooper pairing condensate and chiral condensate at finite baryon density. 

By reformulating the IT model [|10| using the Bogoliubov-Valatin trans- 
formation, we show that the mixed phase, in which both Am and are 
non-zero, is energetically favored. By symmetry we argue that the ISB phe- 
nomena may not be a relevant notion in the IT model. 

^If we consider density dependence of mp and A^, the fermion number density may 
not be zero but is still expected to be small. 
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0.3 



0.2 - 



0.1 



0.2 0.4 0.6 0.8 1 

Figure 2: Fermion number density as a function of chemical potential. Here 
p = and /2 = . In the regime /i < /i^ the system sits in the chiral symmetry 
breaking phase(J/) and in the mixed phase phase(/) in the case of /i > /ic- 



We calculated the effective potentials V{4> = 0, Vo,M) and V{M = 
0,Vq,(P) which are functions of the fermion-antifermion (or chiral), ISB (or 
mean-field) and Cooper-pairing order parameters (/>, Vq and M, respectively. 
At zero chemical potential, we find that the global minimum of the effective 
potential is given by Vq = and M / or Vq = and cp ^ 0. Since Mq 
should be zero [15| or very small, it is reasonable to conclude that the system 
sits in the chiral symmetry breaking phase. 

At finite density, taking A/j = vr and a = — O.lvr, we find that in the 
regime fj, < Hc, we have the chiral symmetry breaking phase (II) with almost 
zero fermion number density and in the case of ^ > the system sits in the 
mixed phase with nonzero "vacuum" (or rather ground-state) expectation 
values of (^^V*) (V'V'); ^-C; ^ 7^ and M 0. We observe an ISB-like 
behavior in the fermion number density as in fig. ^ but we fail to observe 
the exclusive competition obtained in Q. As we can see in (^) and also 
in ref.[15|, M (~ (V'V')) is independent of the chemical potential. We are 
unable to say whether or not these are an artifact of the simplified models 
not present in effective theories of QCD. A similar (uncertain) situation 
applies also to analysis made in QCD at weak coupling |[T7[| . 

Our analysis depends on the value of the renomalization point which is 
arbitrary in general. The possible resolution to this arbitrariness is to do 
a renormalization group analysis at finite density discussed in [16|. Such a 
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procedure will replace the renormalization scale (up to a constant) by the 
chemical potential [|l^]. We leave this exercise to a future publication. 
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Appendix 

In this appendix, we compare the grand canonical ensemble used in this 
work and in ref.|jlO|, |ll| (case B) with that widely used in BCS theory, 
for example see (case A). Here (a, a*) stand for the annihilation and 
creation operators for the "bare" fermions and {b,b*) for the quasiparticles 
as in section 2. 

Ab + 0, = 

Here we shall take the mean-field gap to be zero, i.e., A^f = for 
simplicity. Since the BCS ground state does not preserve the particle number 
of the ground state, a condition is imposed on the average; 

{^BCs\Nop\^BCs) = N (72) 

where Nop = J2k'^*k^k, the number operator. As described in ref.jl^ in 
detail, this relaxation of the condition of N in the BCS and Bogoliubov 
theory must be distinguished from the use of grand canonical ensemble in 
statistical mechanics. In grand canonical ensemble of statistical mechanics, 
we deal with an ensemble of systems, with a distribution of particle numbers 
N with the systems with N particles being weighted by a factor where z 
called fugacity is defined by z = e^^. However, each separate system within 
an assembly has its own definite number of particles N. 

In the case of A, we add a term —^a(^*ci to the Hamiltonian to specify 
the mean value of particle number. Normally in this case, we adjust the 
chemical potential /i to satisfy {Na) = Na- Here the chemical potential is 
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a variable conjugate to Na = J2k (^*k^k- After specifying the mean number of 
Na, we diagonalize the Hamiltonian to get H(b*b) describing quasi-particles. 
We can specify a constant number (in the sense of the average given above) 
to a system characterized by this diagonahzed Hamiltonian H(b*b). There- 
fore this system becomes effectively canonical ensemble of the statistical 
system with definite number but we are essentially using the grand canon- 
ical ensemble. Then the probability of quasi-particle thermal excitation is 
given at a given energy e by (6*6) = (1 -|- e^v^^^~'^°^^+^)~^. 

In the case of B, we do not constrain the mean number in terms of Na = 
J2k ^*k^k- Since our diagonahzed Hamiltonian H(b*b) has the same structure 
with that of non-interacting harmonic oscillator and preserves the particle 
number, we can define a simultaneous eigenstate of quasi-particle number 
and energy operators 0. Note that the quasiparticle operator b kills, by 
definition, the BCS ground state 6|^'bcs) = and the lowest quasiparticle 
excitation state is given by b*\^BCs)- Following an elementary procedure in 
statistical mechanics, we have {b*b) = (e^(V^^+^s~'^*) + with being 
the "chemical potential" for quasiparticles. This formalism B is essential if 
we want to consider both > and < on the same footing. If we 
use the formalism A, the gap equation becomes 

/ /? J(e - fia)^ + A% 

^(e-/.J2 + A| = ABtanh-V: ^ 1, (73) 

Taking /? ^ 00, we have 

^{e-fiar + Al = XB (74) 

showing that Ab < does not enter into the consideration. Note that Xb > 
corresponds to the usual BCS attraction. 

Since the chemical potential fib essentially controls the number of quasi- 
particle excitations while the chemical potential fia fixes the mean number 
of the ground state, one might say that fib has nothing to do with fia- But 
one should keep it in mind that b* is a linear combination of the a and 
a* original fermion operators. Thus b* can create and annihilate fermions 
{a, a*) from the ground state. This implies that controlling the number of 
quasi-particles must have something to do with that of the "bare" fermion 
(a*, a). Whether we fix the number of ground state(^) or quasi-particle 

®This is not true in the case of BCS and Bogoliubov theory. 
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excitation(i3), they should give the same answer for {a*a). Note that (a*, a) 
and {b*,b) are connected by the canonical transformation: 

= + T^B^} with ii; = - f^a? + A| (75) 

where e = p^. 



Ab = 0, Am^O 

In this case the probability of thermal excitation of "bare" fermion is 
given by (for e + Am > fJ-a) 

(a*(p)a(p)) = a/'^, A r with B 

with A 



1 _|_ g/3(e+AM-Ma) 

= {h*{p)h{p)) (76) 

showing = ^^b■ Note that in the case = 0, / 0, the KMS-state 
is defined for the operator (a,a*)|n 
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